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its  practical  implementation  from  all  the  methods  mentioned.  This  radical 
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6c  terminated,  having  found  an  invariant  aubapace.  I  g_2.  The  Lon 


ant*  of  the  Lanczot  algorithm  in  order  to  determine  among  »ev*ral  poasible 


Denote  by  X,(IT)  the  emelle.t  end  by  X,(»)  the  t*r*eet  eigenvalue  of  the 

matrut  W. 


Then  considering  the  k-th  step  in  the  net  way  as  the  first  step,  one  finds 
that  3k  end  can  b®  computed  by  the  recurrence 


quence  rather  than  a  cause  of  the  loss  of  orthogonality,  and  will  not  concern 


Ufl  for  the  moment  norm  of  this  perturbation  can  be  bounded  by 


turbation  of  Ap  »  NfANf. the  orthogonal  projection  of  A  onto  span(Qj). 


orthogonality  at  a  level  of  t.  then  the  Lanczos  algorithm  actually  computes  a 


*.6  SKniorthoceoalilaUea 3tr*te*te*.  j)  The  number*  i«l.  ■  ■  J-l  are  chosen  *uch  that 

A*  *oon  u  the  leeel  ot  orthogonality  deteriorate*  »o  much  that  l«V»tl  *  |«,„,|  *  u»  (2.1.2) 

IfAid*  I  >  «*  •  tor  some  k  <  j.  the  nice  result  T^NjANj  doe*  not  (2.1.2) 


iiiiona  are  aatial 


Uu<  method  one  h«i  to  update  en  n 
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orthogonality  arc  plotted  In  figure  3.1. 


-  Estimated  level  of  orthogonality.  «  =  1000.0 


reach  the  critical  Vt  region.  Even  the  curve  with  the  largest  overestimate 
reaches  the  Vc  threshold  only  three  steps  to  early  at  step  6Q.  In  spite  of  the 
dependence  on  the  random  terms.  (3.1.2)  therefore  appears  to  produce  a 


(31.2)  signals  at  about  the  right  Lanczos  step  that  the  V* -level  has 


•Tftnd.dev.  1  no.  of  itc nil 


This  information  about  |/j  I  i*  nr y  helpful  for  the  study  of  f//,  -  q£/i.  ,  8m 

=  - VtVj-i  *  Oj  .  (3.1.8) 

Since  fj  and  fa  ere  unit  vectori,  If)  and  I/*  we  expect  Pl" 

where  if  i>  the  turn  of  all  remaining  terms.  We  therefore  computed  the  term 
that  (f//e  -  lift)/  Thi*  expectation  proved  to  be  correct.  , 

Th-  reaulte  of  the  text  rune  with  the  tame  text  matrices  as  before  are  sum-  ^'filin''  The  "tr*  n  *U  l»tr«lue«d.  since  we  are  trying  to  estimate  the 


Since  the  two  vector*  were  already  semiorthof  onal  before  the  reorthof  onali-  5,2.  Tha  Behavior  of  the  Coen  paled  Uml  of  Orthogonality- 

»Uon.  Ih.  nlue.  in  T.bl.  3.4  »re  quit.  »>U.i>d  there  i»  no  .Iron*  dep.n-  ln  ieelion  3. ,  „  ai,cue,.d  bow  en  eetimeU  for  the  lew.)  of  ortboqoMlity 

dence  on  n.  Beeed  on  tbie  reeult  -e  .re  cbooemg  u,.,.«(".15)>  efter  e  CU1  be  computed  with  formal.  (3.1.2)  end  .ppn.pri.tely  cboMn  r*ndom 

reortbogoraliution  hu  been  performed.  number,  for  dp,  (3.1.3).  ft  (3.1.7).  end  u*  »fter  reorthOfW»lU»tion  (3.1.8) 
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the  tire  of  f/’.rft-  If  bow  ever  two  reortbogonalitation#  are  performed  In  a 
row.  then  formula  (3.1.1)  field# 


figure  3.8  »bow»  the  typical  pattern  in  the  loaa  of  orthogonality  Ueually  only 


some  neighboring  g/+i9*  hove  grown  to  about  the  V?  level,  whereas  most  against  therefore  seems  to  group  them  into  "batches".  One  batch  contains  ail 

other  inner  products  remain  quite  small  lo  order  to  maintain  semiortho-  the  offending  Lancsos  rectors,  i.e„  all  the  9*  with  I9j%t9* !  >  *nd  in  addi- 

gonality  it  is  therefore  only  necessary  to  octhogonallie  against  selected  tion  to  kbat  m  certain  number  of  neighboring  vectors.  The  nest  obvious  ques- 

Lanczos  vectors  In  the  example  given  in  the  table  it  could  be  the  first  ten.  tion  ia  then:  bow  many  neighboring  vectors  should  ba  included  in  those 


I 


best  strategy  for  chosing  Lancrot  vectors  to  reorthogonalize 


border  the  batch.  because  we  do  not  orthogonalize  against  them  at  the 


0+iy-*t  atop  anyway.  Therefore  it  U  poaaibla  to  compute  already  at  the  J-th  |  J.B.  Compartaoo  with  Select!*. 


A  «  ding  (100, 48.5, 47.5...  .-47.5. -48.5.-49.5)  and  the  riiht  hand  »ide  If  we  consider  on  the  other  hand  elf envalue  computations,  then  the  orl- 

6  =  (100.48.5.47.5,  .-4 7. 5, -♦8.5. -49. 5).  Solving  Ax  *  b  we  encounter  the  glnal  problem  can  be  changed  by  shifting  and  inverting  In  such  a  way  that 

special  situation  that  the  solution  vector  has  equal  components  in  direction  the  values  and  vectors  which  converge  first  are  indeed  the  desired  ones, 

of  all  eigenvectors,  and  that  the  eigenvalues  of  A  are  evenly  distributed  (with  jn  cose  (compare  example  3)  SO  fares  much  better,  since  the  loss  of 


] 


[8]  G.  Golub,  R.  Underwood,  and  J.H.  WUkinton.  The  Lanczo*  Algorithm  for  (lBl  c-  Uneio*.  Solution  of  Syeteme  of  linear  Equation,  by  Minimized  Here- 

the  Symmetric  Am  e  \Bx  Problem.  Tech.  Rep.  STAN -CS- 72- 720.  Comp.  tion *• 1  *»*•  Bur-  Standard.  4®.  33-53.  1P52. 

Set  Dep..  Stanford  Unleenity,  1872.  [18]  D.C.  Luenberger,  Hyperbolic  Poire  In  tbe  Method  of  Connate  Crediente. 

[8]  l.  Greer.  Analyaaa  of  the  Lancxoa  Algorithm  and  of  the  Approximation  S1A1* 1  AppL  17-  1253-1287.  1888. 


This  report  was  dona  with  support  from  the  Cantor 
for  Part  and  Applied  Mathematics.  Any  conclusions 
or  opinions  expressed  In  this  report  represent  solely 
those  of  the  author(s)  and  not  necessarily  those  of 
The  Regents  of  the  University  of  California,  the 
Center  for  Pure  and  Applied  Mathematics  or  the 
Department  of  Mathematics. 


